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An integrated approach for entry trajectory design, guidance, and simulation 
is proposed. The key ingredients for this approach are an on-line 3 degree-of- 
freedom entry trajectory planning algorithm and the entry guidance algorithm that 
generates the guidance gains automatically. When fully developed, such a tool 
could enable end-to-end entry mission design and simulations in 3DOF and 6DOF 
mode from de-orbit burn to the TAEM interface and beyond, all in one key stroke. 

Some preliminary examples of such a capability are presented in this paper that 
demonstrate the potential of this type of integrated environment. 


I. Introduction 

The traditional approach for entry mission plan- 
ning and trajectory design is based on numerical 
optimization methods. The process is laborious 
and requires familiarity with the optimization soft- 
ware and experience with the optimization method 
because man-in-the-loop adjustments are inevitably 
needed. Consequently, trajectory analysis and de- 
sign are typically restricted to baseline missions when 
only a limited amount of time is available. This pro- 
cess is also separated from the simulations needed 
to verify that the guidance and control systems can 
actually fly the designed trajectory. These verifica- 
tions and validations happen only downstream when 
the reference trajectory data has been processed and 
translated into appropriate guidance system I-loads. 
A new approach based on a recently developed al- 
gorithm for entry trajectory planning and guidance 
has the potential to significantly shorten this time- 
consuming process and allow quick turn-around in 
mission planning and trajectory design. 1 The algo- 
rithm is originally developed with the objective of 
on-board applications. But the efficiency and versa- 
tility of the algorithm also make it ideal as the core of 

* Professor, Department of Aerospace Engineering, Asso- 
ciate Fellow AIAA. Email: plu@iastate.edu 

t Advanced Space Programs, Associate Fellow AIAA 
Copyright (c) 2004 by the American Institute of Aeronautics 
and Astronautics, Inc. All rights reserved. 


a tool that enables rapid entry trajectory design and 
mission planning with 3DOF or 6DOF simulations 
all at once. This algorithm has two parts. The first 
is a 3DOF entry trajectory planner that generates 
within 2-3 seconds on a desktop computer a feasible 
entry trajectory based on the vehicle model, current 
state, targeted final conditions and trajectory path 
constraints. The second is a self-tuned entry guid- 
ance law that automatically computes the guidance 
gains based on the reference trajectory generated by 
the trajectory planner. Therefore no post-processing 
and generation of guidance I-loads are required af- 
ter the reference trajectory is found. The algorithm 
is readily implementable in a simulation environ- 
ment where vehicle dynamics and other subsystems 
such as control, navigation and avionics are modeled. 
In this environment end-to-end simulations of entry 
flight can be performed with just one key stroke. 
Therefore many different possible mission scenarios 
involving different landing sites with various down- 
range and crossrange conditions can be investigated 
in trade studies in a relatively short period of time. 
The effects of vehicle parameters such as lift-to-drag 
ratio, angle of attack profile, and mass as well as 
trajectory constraints in peak heating rate or tem- 
perature, deceleration and dynamic pressure can be 
rapidly assessed. The information gathered in such 
studies may be passed back to the vehicle design 
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team if necessary in the case of development of a 
new entry vehicle. 

This paper will demonstrate some of the above 
mentioned capabilities of this technique. The mis- 
sion scenarios are for an entry vehicle returning from 
the International Space Station orbit. The vehicle 
has a medium lift-to-drag ratio, and is to land at the 
Kennedy Space Center. The cases studies involve 
different entry conditions resulted from different de- 
orbit opportunities, several nominal angle of attack 
profiles, and parametric variations in vehicle aerody- 
namic properties. 


II. Entry Trajectory Planning Algorithm 

The reference trajectory for the Space Shuttle is 
represented by a drag acceleration profile planned on 
the ground . 2 By tracking this drag profile onboard, 
the longitudinal parameters of the entry trajectory 
(altitude, velocity, and range) are determined. The 
lateral motion is controlled by a bank reversal logic 
based on heading angle offset to the landing site. 
Different approaches have continued to be be inves- 
tigated since the first flight of the Shuttle . 1 ’ 3-8 The 
more recent efforts have focused on developing the 
capability of planning the reference trajectory on- 
board. 

At the core of this fast approach for entry mis- 
sion analysis and evaluation is a fast and reliable 
trajectory planning algorithm. Given initial condi- 
tions, landing site coordinates, terminal constraints, 
and trajectory path constraints, this algorithm will 
generate within 2-3 seconds a 3DOF feasible en- 
try trajectory. This algorithm, implemented in a 
simulation environment, allow seamless integration 
of trajectory planning, guidance, control, and other 
subsystems in the simulations. This trajectory plan- 
ning algorithm is described in detail in Ref. 1. For 
completeness and the convenience of the reader, this 
algorithm is briefly reviewed in this section. 
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where r is the radial distance from the center of 
the Earth to the vehicle, normalized by the ra- 
dius of the Earth Ro = 6378 km. The longitude 
and latitude are 0 and <p, respectively. The Earth- 
relative velocity V is normalized by V c = y/goRo 
with g 0 — 9.81 m/sec 2 . The terms D and L 
are the aerodynamic drag and lift accelerations in 
g' s, i. e., D = p(V c V) 2 S re fC D / (2mgo) and L - 
P(V C V) 2 £re/C f z,/( 2 mpo)j where p is the atmospheric 
density, S re f the reference area of the vehicle, and 
m the mass of the vehicle. Note that D and L 
are also functions of a, the angle of attack, through 
the dependence of the drag and lift coefficients C D 
and Cl on a. The flight path angle of the Earth- 
relative velocity is 7 and a the bank angle. The 
velocity azimuth angle xp is measured from the North 
in a clockwise direction. The differentiation is with 
respect to the dimensionless time r — t/y/Ro/go - 
Final ly, Q. is the Earth self-rotation rate normalized 

by Vgo/Ro ■ 

Typical inequality entry trajectory constraints in- 
clude: 


Q < Qmax 

(7) 

|Lcosa + Dsina| < n Zmax 

( 8 ) 

Q ^ Qmax 

(9) 


2.1 Entry Dynamics and Trajectory Con- 
straints 


i — V 2 — L cos a eq < 0 
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The 3DOF point-mass dynamics of the entry vehi- 
cle over a spherical rotating Earth are described by 
following dimensionless equations of motion 9 

r = Vsin 7 (1) 
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where Eq. (7) is a constraint on the heating rate 
at a specified stagnation point on the surface of 
the vehicle, with Q = ky/pV 31 5 for a constant k . 
When applicable, this constraint may be replaced by 
a altitude-versus-velocity curve that represents the 
temperature limits at various body points. The con- 
straint Eq. ( 8 ) is on the aerodynamic load (in g's) 
in the body-normal direction. Depending on the ve- 
hicle configuration and mission, this constraint may 
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be replaced by the total load constraint: 


V L 2 + D 2 < n max (11) 

Constraint Eq. (9) is on the dynamic pressure 
with q — p(V c V) 2 /2. The parameters Q max , 
n Zmax( n max)i and qmax are vehicle-dependent con- 
stants. The constraints in Eqs. (7)-(9) are the 
“hard” constraints, meaning that they should be 
enforced strictly. When held in equality, the last 
constraint Eq.(10) is called the equilibrium glide con- 
dition (EGC) at a ~ a eq 7 where oeq is a specified 
constant. The EGC is obtained by omitting the 
Earth self-rotation terms and setting 7 = 0 and 
7 = 0 in Eq. (5). This constraint serves to reduce 
the phugoid oscillations in altitudes along the en- 
try trajectory, and preserve bank angle margin so 
that cr > a eq during most of the flight. The lat- 
ter allows sufficient trajectory control to account for 
dispersions. The constraint Eq. (10) is a soft con- 
straint in the sense that its observance need not be 
absolutely strict, particularly in the final portion of 
the entry trajectory. 

The entry trajectory terminates at a condition 
traditionally referred to as the Terminal Area En- 
ergy Management (TAEM) interface. This condition 
can be for instance a given value of velocity /Mach 
number. At the TAEM interface, the entry trajec- 
tory must have correct conditions to ensure success- 
ful TAEM flight and landing that will follow entry. 
These conditions form the terminal constraints on 
the entry trajectory in terms of: 

*7 = r T AEM , Vf — Vtaem , $/ = staem ( 12 ) 

where s/ is the value of range- to-go s t0ff0 , defined 
to be the range from the vehicle to the tangency 
of the heading alignment cone (HAC) near an end 
of the runway for horizontal landing. The coordi- 
nates of the center of the HAC are known. In the 
case of other forms of landing and recovery (e.g., 
via parachute), HAC will represent the landing site 
although the term “HAC” is still used in this pa- 
per in a generic sense. The TAEM altitude vtaem , 
velocity Vtaem and range- to-H AC staem are all 
specified for a given vehicle. In addition, the Earth- 
relative velocity vector at the TAEM interface should 
be pointing nearly to the HAC tangency. Let A tj)f 
be the difference between the velocity azimuth angle 
and the line-of-sight angle from the vehicle to the 
HAC at the TAEM interface. This condition is then 
expressed by 

(13) 


This condition stipulates that the final velocity vec- 
tor should be directed at the HAC within a given 
tolerance A iPtaem- 

The bank angle magnitude at the TAEM interface 
often times is also another parameter that is desired 
to be constrained. Too large a a f could result in large 
transient response for TAEM guidance and control, 
and even make recovery unachievable for a less ma- 
neuverable entry vehicle. Thus, the constraint 

Wf \ < ^taem (14) 

for a given otaem > 0 may also be imposed. 

2.2 Trajectory Planning Problem and Algo- 
rithm Outline 

The entry reference trajectory generation prob- 
lem is defined as follows: given the entry conditions, 
the path constraints and terminal conditions at the 
TAEM interface as in the preceding section, find 
the state history of x = {r 0 <j> V J ip} T and the 
corresponding trajectory control u = {a a} T such 
that:the state and control profiles x(£) and u (t) sat- 
isfy the 3DOF equations of motion Eqs. (1-6), all 
the required TAEM interface conditions, and the im- 
posed trajectory path constraints. A pair [x(£) u (£)] 
that meets the above conditions is called a feasible 
trajectory. We are concerned with only lifting entry 
flight, i.e., L/D ^ 0. Furthermore, we will assume 
that a nominal angle of attack profile is also specified 
for a given vehicle as it is in many cases. Such an 
a- profile usually is determined on the basis of flight 
trim conditions and thermal considerations. 

The trajectory planning algorithm developed in 
Ref. 1 tackles the problem in two major sequen- 
tial steps: first determining the longitudinal profiles 
of the trajectory, and then completing the lateral 
profiles. In each step there is one single-parameter 
iteration involved. All the trajectory constraints ex- 
cept for Eq. (13) are enforced by the longitudinal 
profiles. The condition (13) is met by finding the 
point for the last bank reversal maneuver, and the 
equations of motion (1-6) are numerically integrated 
in the second step. 

In the first step, the whole entry trajectory is 
divided into three phases. Each phase has its distinc- 
tive physical characteristics and the algorithm uses 
different approach in different phases: 

1. Initial descent phase: the flight in this phase is a 
“controlled fall” in which the vehicle uses a con- 
stant bank angle, and the trajectory takes the 
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vehicle from the entry interface at about 120 km 
in altitude to an altitude at about 80 km. This 
is where the dynamic pressure has reached a suf- 
ficient level for aerodynamic lift to become influ- 
ential in shaping the trajectory. The magnitude 
of the appropriate constant bank angle is deter- 
mined by increasing the bank angle from zero at 
a fixed incremental and numerically integrating 
the equations of motion. The incremental pro- 
cess is stopped when the condition at the end of 
the initial descent satisfies a smooth transition 
into the next phase of quasi-equilibrium glide. 

2. Quasi-equilibrium glide (QEG) phase: The 

QEG phase covers the majority of the entry tra- 
jectory where all the inequality path constraints 
must be observed and the range achieved must 
be correct for the vehicle to reach the landing 
site. In this phase the Sight path angle 7 is small 
and varies relatively slowly. Setting cos 7 = 1 
and 7 = 0 in Eq. (5) and ignoring Earth rota- 
tion gives 

Lcoso + (V 2 — —)— = 0 (15) 

r r 

This relationship is called the quasi-equilibrium 
glide condition (QEGC). It is valid in lifting 
entry flight from an altitude of about 75-80 
km to certain supersonic velocity which is de- 
pendent on the lifting capability of the vehicle. 
The QEGC is used in the trajectory planning 
algorithm as a motion integral. It allows conve- 
niently the bank angle to be tied directly to the 
longitudinal trajectory state variables velocity 
and altitude. With the aid of the QEGC, the 
satisfaction of the requirements on range, veloc- 
ity and altitude and observation of the inequal- 
ity trajectory constraints are all accomplished 
by the iterative determination of a bank angle 
parameter (typically the value of bank angle at 
the midway point along the QEG phase). 

3. Pre-TAEM phase : In this phase the velocity has 

reduced to a level close to the TAEM velocity, 
and the flight path angle as well as the altitude 

starts to decrease more rapidly. The QEGC is 

no longer valid. The algorithm numerically in- 
tegrates the equations of motion backward from 
the given TAEM conditions to a selected pre- 
TAEM velocity (usually around Mach 5) where 
the QEGC phase connects with the pre-TAEM 
phase. The bank angle is defined by flying an 
altitude profile as a polynomial in velocity. No 
iteration is involved in this phase. In terms 
of execution sequence, the pre-TAEM phase is 


planned before the QEG phase so that the pre- 
TAEM results establish well defined conditions 
at the end of the reverse integration where the 
QEG phase ends. These conditions are then 
used in determining the bank angle profile in 
the QEG phase. 

In the second major step, the longitudinal profiles 
(range, altitude, and flight path angle versus energy) 
found in the first step are tracked with linearized 
time- varying feedback control laws in the magnitude 
of the bank angle and angle of attack when angle 
of attack modulation is allowed. The sign of the 
bank angle is initially such that the heading error 
to the HAC is reduced. If the sign of the bank 
angle remains unchanged, the trajectory would over- 
correct itself and miss the landing site. Therefore 
the sign is reversed at an appropriate point so that 
when the vehicle reaches the TAEM interface, the 
condition Eq. (13) is satisfied. The determination 
of the correct bank reversal point is the second one- 
parameter iteration in the algorithm, and involves 
numerically integrating the full equations of motion 
with the bank angle and angle of attack from the 
tracking laws. Once this second search is finished, 
a complete 3DOF feasible entry trajectory is found, 
and ready to be used as the reference trajectory for 
entry guidance. 

III. Entry Guidance Algorithm 

The entry guidance law used in this paper is based 
on receding-horizon control (RHC) approach. Con- 
sider the linearized dynamics of the system Eqs. 
(1-6) about the reference trajectory 

<5x = A(t)<5x + B(t)6\ 1 (16) 

where <5x € R 3 * * 6 is the state dispersion, Su — (So 6a) T 
the control dispersion, and A(-) : R — ► R 6x6 and 
B(-) : R — > R 6x2 are the Jacobians of the system 
Eqs. (1-6) with respect to the state and control 
vector, respectively. These two matrices are time- 
varying because the reference trajectory is time- 
varying. The regulation/stabilization problem seeks 
to find a control law Su = K(t)Si c to drive <5x — > 0 for 
any initial nonzero condition <5x(to)- The receding- 
horizon control problem is defined to be an optimal 
control problem in which the performance index at 
any t 

[t+T 

J — J [Sx T (r)QSx(r) + Su T (r)RSu(T)\dT (17) 

is minimized for some Q > 0 and R > 0, subject 
to system dynamics (16), initial condition 5x(t), and 



the terminal constraint 

Jx(< + T) = 0 (18) 

where 0 < T < oo is a constant that defines the con- 
trol horizon of the problem. The idea is to solve this 
optimal control problem in the finite moving horizon 
[t, tAT] with the current state x(t) as the initial con- 
dition. Let Suoptm (*) be the optimal control obtained 
in [t, t A T]. Only the first data Su optm (t) is used to 
be the current applied control Su(t) ~ £u op * m (i). 
The rest of Su op tm{’) is discarded. The process is 
then repeated for the next t. Because the control 
<5u(-) applied at any t is dependent on the condition 
<5x(£), the receding-horizon control strategy produces 
a feedback control. 

A closed-form approximation to the receding- 
horizon control law has been developed. 10 Consider 
the above receding-horizon problem in the interval 
[t, t 4- T\. Divide this interval into TV subintervals of 
equal length h — T/N for some integer TV > 3. With 
the current state <5x(i) known, a first-order predic- 
tion of 6x(t -I- h) as a function of Su(t) is given by a 
Taylor series expansion at t 

Sx(t 4- h) « 6x(t) + /i[A(£)£x(t) + B(t)Su(t)] 

= (I A hA)8x A hBSu (19) 

Denote A k = A(t + kh), Bk = B(t A kti), Sxk — 
8x(t A kh) and Su k — 8u(t A kh ), k = 1, TV. Then 
another first-order Taylor series expansion at t Ah, 
together with Eq. (19), gives 

Sx(t A 2h) « 8x i A h[Ai 5xi A i?i$ui] 

« (I A hA x )(I A hA)5x 
A h(I A hA\)B6\x A hBrfm (20) 

Continuing this process, we have 
k - 1 

Sxk Rs Afc<5x + ^2 k = l,...N (21) 

t=0 

where 

A fc = (I A hA k „ i)A fc _x, with A 0 = / (22) 
Gk,i = (7 A hAk—^Gfc—i'i, i = 0, 1, ..., k — 2, 
and i — hB k —i (23) 

The subscript 0 in above expressions denotes the val- 
ues at t. Let L 0 — 5x r (t)Q^x(t)-hJu r (t)Mu(t) and 
L k = Sx k Q5x k A Su^R8u k , k = 1, ..., TV. The in- 
tegral in Eq. (17) is approximated by the standard 
trapezoidal formula for integrals: 

J ~ /i(0.5Lo A L\ A ... + L,/v_i + 0.5 Ljy) (24) 


Define an (2TV)-dimensional vector v = 
co/{£u(£), <5ui, ..., 6u n _i}. If the 8x k s in Eq. (24) 
are replaced by the prediction (21), the performance 
index is thus approximated by a quadratic function 
of v: 

J = ^v T H(t,N,h)v + 6x T S(t,N,h)v 

A q{8x , t, TV , h, Sun) (25) 

where H € R 2Nx2N is positive definite for any t > to , 
integer TV and h > 0, S G R 6x2N and q is quadratic 
in <5x and <5u/v. These terms are obtained directly 
by rearranging the expression of J in Eq. (24) once 
8x k s are replaced by Eq. (21). Examples of the 
expressions of H and S for several values of TV can 
be found in Ref. 10. The constraint Eq. (18) can 
be approximated by setting Sx n = 0 from Eq. (21), 
which can be rewritten in a compact form 

M r (t, TV, h)v — - A n 6x (26) 

where 

M t = [Gw, o ... Gjv,n-i] € R 6x2N (27) 

The minimization of J in Eq. (25) with respect to 
v subject to constraint (26) constitutes a quadratic 
programming (QP) problem. When M T has full 
rank, the above QP problem has a unique solution 

v = -{[if -1 - H- 1 M(M t H~ 1 M)~ 1 M t (28) 

x zr 1 ] S T + [i/" 1 M(M 7 ’ff- 1 M) -1 ] Aat}<5x 

Define an 2 x 2TV matrix 

hN = {^2x2, 0, ....0} (29) 

A closed-form, linear, time-varying feedback control 
law for <5u(t), denoted by Su(t;N,h) to signify its 
dependence on time and the values of TV and h , is 
then obtained from the first two equations in Eq. (28) 

Su(t; TV, h) = I 2 nv = K(t, TV, h)Sx(t) (30) 

The guidance law Eq. (30) regulates the trajectory 
in both longitudinal and lateral directions to follow 
the reference trajectory. All the feedback gains are 
generated automatically from the vehicle model and 
reference trajectory as outlined above. Therefore 
there is no need for a separate bank reversal logic 
in the entry guidance algorithm as it is in the Shut- 
tle entry guidance. This is particularly convenient 
during the development and analysis phase when the 
vehicle design is evolving, and different mission sce- 
narios are examined. No extensive guidance logic 



re-design and adjustments are required for each new 
vehicle design update and changing mission profile. 
However, this approach has one possible drawback 
in the operational phase when significant aerody- 
namic modeling uncertainty exists. In the presence 
of large aerodynamic modeling mismatch, it might 
become difficult for the vehicle to follow closely both 
the longitudinal and lateral profiles of the reference 
trajectory which is designed on the basis of the nom- 
inal aerodynamic model. In such a case a guidance 
law may be more robust that only flies the reference 
longitudinal profiles, and controls the heading of the 
vehicle by a simple bank reversal logic. 6 

IV. Applications 

To demonstrate the application of the the above 
described approach, entry missions from the Inter- 
national Space Station (ISS) orbit for landing at the 
Kennedy Space Center (KSC) are considered. The 
entry vehicle has a medium lift-to-drag ratio. The 
TAEM interface is defined at a relative velocity of 
Mach 3. The following TAEM conditions are im- 
posed: 

altitude - 24.384 km (80,000 ft) (31) 
range-to-HAC = 55.56 km (30 nm) (32) 
heading-to-HAC offset < 5 deg (33) 

In addition, the trajectory is further constrained by 

Q < 1134893 W/m 2 (100 BTU/ft 2 -sec) (34) 
y/l , 2 + D 2 < 2.5 g (35) 
q < 14364 N/m 2 (300 psf) (36) 

The simulation environment is a tool developed 
at the NASA Marshall Space Flight Center, called 
Marshall Aerospace VEhicle Representation In C II 
(MAVERIC II). This is a modularized high-fidelity 
simulation software that allows the user to use ei- 
ther built-in sub-system models or replace them 
with different vehicle model, sub-systems and algo- 
rithms. The mean and dispersed wind data from 
the Global Reference Atmospheric Model (GRAM) 11 
is available in MAVERIC II for nominal and dis- 
persed simulations. The entry guidance algorithms 
described in the preceding sections are implemented 
in MAVERIC II for the vehicle used in this paper. 
All the simulations are performed in MAVERIC II. 
At this point the control system design is not yet im- 
plemented for this vehicle, thus all the simulations 
are 3DOF simulations. 


4.1 Entry Flight Assessment for Multiple De- 
Orbit Options 

The objective of this type of analysis is to assess 
the capability of the vehicle to accommodate a range 
of de-orbit opportunities and still ensure a successful 
entry flight toward the primary landing site. The or- 
bit from which the entry vehicle returns is assumed 
to be an orbit with the same altitude and inclination 
as those of the ISS orbit. Six de-orbit opportunities 
are examined here that result in different entry con- 
ditions at the entry interface defined at the altitude 
of 121.92 km (400,000 ft). These cases correspond to 
tests 13 to 17 in Ref. 12. The entry conditions are 
summarized in Table 1. They are separated into two 
groups, Mission A1 - Mission A3, and Mission B1 - 
Mission B3. The defining difference between the two 
groups is in their entry flight path angles: group A 
has a shallower entry night path angle thus longer 
downranges, whereas group B has a steeper entry 
flight path angle and shorter downranges as a result. 
In each group, the first mission (A1 or Bl) corre- 
sponds to de-orbit at a best opportunity for landing 
at the KSC. Therefore the crossrange is minimum. 
The other two missions in each group are for the 
cases when the de-orbit maneuver takes place one 
orbit earlier and one orbit later than the best oppor- 
tunity, respectively. The consequence is a large right 
crossrange or a large left crossrange at the entry in- 
terface. The negative sign for crossrange in Table 1 
indicates a left crossrange. 

The integrated approach described in this paper 
enables an end-to-end simulation for each of the pos- 
sible de-orbit opportunity without the need to wait 
for off-line generation of a reference entry trajectory 
for a particular set of entry conditions. Instead, the 
reference entry trajectory is planned in the first en- 
try guidance cycle. If the entry conditions are such 
that a feasible reference trajectory cannot be found 
for the vehicle, the trajectory planner sends back a 
flag to so indicate and provide a possible reason why 
a trajectory does not exist. The de-orbit maneuver 
then should be adjusted accordingly. 

All the six set of entry conditions are feasible for 
this entry vehicle. Figure 1 shows the ground tracks 
of the 6 missions. Figures 2 and 3 illustrate the com- 
parisons of altitude- vs- velocity, bank angle and angle 
of attack histories of the reference trajectory and the 
simulated trajectory for Mission Al. The closed- 
loop guided trajectory follows the reference trajec- 
tory quite well in most of flight. The differences in 
initial portion of the trajectory can be attributed to 
the differences between the high-fidelity simulation 



environment (including winds) and the relatively less 
extensive modeling used in the trajectory planning 
algorithm for efficiency. The same comparisons for 
Mission A2 and A3 are included in Figs. 4 and 5 
where the similar observations can be made. Figure 
5 exhibits close symmetry in bank angle histories 
for Mission A2 and A3, which is to be expected, 
given that the entry conditions for these missions 
are nearly symmetric. This result also demonstrates 
the desired predictability of the outcome of the tra- 
jectory planning algorithm. 

Table 2 summarizes the TAEM condition precision 
of all the 6 missions from MAVERIC II simulations. 
The precision is easily at the best possible level as 
far as the performance of any type of entry guidance 
is concerned. 

4.2 Flying Different a-Profiles 

It is well known that the angle-of- attack profile 
flown by the entry vehicle affects considerably the 
footprint of the entry trajectories. Perhaps some- 
what less well known is the fact that different a- 
profiles also have different predictable impact on the 
heat load absorbed by the vehicle. In their classic 
work Allen and Eggers show that the total heat load 
Q can be approximated by 13 



where Cf is the equivalent skin-friction coefficient, 
S w the wetted area, and Vo the entry velocity. Differ- 
ent a-profiles will result in difference values for Co, 
hence different heat load, assuming that the variation 
in Cd does not lead to the variation in C/. Higher 
a will produce higher Cd, and thus lower heat load, 
and vice versa. If the flight times do not differ appre- 
ciably with different a-profiles (as usually the case), 
the heat rate along the trajectory will also follow the 
same trend as the heat load decreases or increases 
when the a-profile is changed. Possible restricting 
factors that may limit the option of flying different 
a-profiles include trim requirement in flight control 
and thermal protection considerations. 

The change of the nominal a-profile will alter the 
entry trajectory, and a new reference trajectory has 
to be planned in a traditional analysis environment. 
The integrated approach presented in this paper 
readily handles the need to assess the effects of flying 
different a-profiles without any additional prepara- 
tion work. The nominal a-profile is an optional input 
to the software, and the simulation is ready to go 
as soon as the a-profile is selected. As a demon- 


stration, a lower a-profile and a higher a-profile are 
used in addition to the nominal a-profile shown in 
Figs. 3 and 5. The difference is a ±5 deg in a at 
hypersonic speeds. The peak heat rate limit for the 
lower a-profile is increased to 1,191,637.7 W/m 2 (105 
BTU/ft 2 -sec), and decreased to 907,914.4 W/m 2 (80 
BTU/ft 2 -sec) for the higher a-profile. All the other 
conditions /constraints remain the same as in Mission 
Al. 

Figure 6 shows the variations of bank angles and 
the 3 a-profiles from the MAVERIC II simulations. 
The a-profile labeled as “medium” is the one used 
in Figs. 1-5 earlier. The trajectories in the velocity- 
altitude space are plotted in Fig. 7. The effects of 
a-profile at hypsersonic speeds are visible in Fig. 7: 
the trajectory is higher with a higher a-profile and 
lower with a lower a-profile. The impact on the heat 
rates is more evident in Fig. 8. Even with the in- 
creased peak heat rate limit, the trajectory with the 
lower a-profile still rides on the heat rate constraint 
boundary in an extended period. On the other hand, 
with a higher a-profile, the trajectory can now eas- 
ily accommodate a 20% reduction in the peak heat 
rate limit when compared to the baseline mission. 
The TAEM condition precision in all these cases are 
comparable to the precision seen in Table 2. 

The change of a profile affects the lift and drag 
of the vehicle, therefore could change the feasibility 
of some missions. For instance, Mission A2 becomes 
infeasible for this vehicle when the lower a-profile 
is flown, because the reduction in lift renders the 
crossrange in Mission A 2 too large for the vehicle. 
On the other hand, all 3 missions B1-B3 are not 
achievable by the vehicle when the lower a-profile is 
used, because the energy dissipates too slowly with 
the shorter downranges (cf. Table 1) and reduced 
drag deceleration. 

4.3 Variations in Cl and Cd 

As the vehicle design evolves, the design modifica- 
tions are bound to bring about changes in the aero- 
dynamic properties of the vehicle. For lifting entry 
any increase in the lift-to-drag ratio will only increase 
the capability of the vehicle to shape its entry trajec- 
tory. It is the reductions in the lift-to-drag ratio that 
will cause the reduction in corssrange capability of 
the vehicle. Therefore the design impact on mission 
needs to be assessed as the design evolution pro- 
gresses. The aerodynamic model is a user-supplied 
module to the entry trajectory planning and entry 
guidance software. Once it is provided, the entry 
guidance software is ready to be used in simulations 



Table 1 Entry conditions at the altitude of 121.92 km (400,000 ft) 


Mission 

V 0 (m/s) 

downrange (km) 

crossrange (km) 

70 (deg) 

A1 

7625.0 

8160 

148 

-1.2492 

A2 

7625.15 

8715 

883 

-1.2493 

A3 

7625.72 

8197 

-778 

-1.2490 

B1 

7622.0 

6519 

59 

-1.4379 

B2 

7621.3 

6554 

809 

-1.4380 

B3 

7622.79 

6589 

-848 

-1.4376 


without further updates. 

With the trajectory for Mission A1 in Section 4.1 
as the baseline, two hypothetic cases are presented 
here to show the application of the proposed ap- 
proach: in one case a uniform 20% reduction in the 
lift coefficient Cl is assumed while Co remains the 
same; the other case has a constant 20% increase 
in Cn with Cl fixed at the baseline level. Note 
that these variations are not treated as modeling er- 
rors, and the trajectory planner and entry guidance 
algorithm are fed with the actual changes. Figure 
9 illustrates the corresponding bank angle histories 
along the two trajectories plus the baseline bank an- 
gle for comparison. In both cases the bank angle is 
smaller compared to baseline, because the L/D ratio 
is decreased. The trajectories in the velocity- altitude 
space are shown in Fig. 10. With the reduction in 
Cx, the trajectory is lower as expected. With the 
increase in Co, the trajectory stays more aloft. The 
effects of the different altitudes at the same veloc- 
ity on heat rate are clear in Fig. 11. Again, with 
the increase in Cp, the heat rate is decreased as dis- 
cussed in Section 4.2. Combinations of simultaneous 
increase in Cp and decrease in Cl have also been 
examined. The qualitative trends are the same as 
seen here. 

When the drag coefficient Cp is increased by 20%, 
this vehicle is still capable of flying all the six mis- 
sions listed in Table 1. But when Cl is reduced by 
20%, the vehicle cannot fly the missions with large 
crossanges anymore (Mission A2, A3, B2, and B3). 

V. Conclusions 

The conventional methodology for entry mission 
design, guidance, simulation and verification follows 
a sequential “divide- and-conquer” approach. Recent 
development of on-line entry trajectory planning al- 
gorithm and entry guidance algorithm with auto- 
matic gains generation allows seamless integration 
of these processes in a single environment where tra- 
jectory design, guidance, control, and simulation are 
all performed with a key stroke. This paper presents 


some preliminary examples of such capability in a 
high-fidelity simulation environment. With the com- 
pletion of development of models and algorithms for 
other sub-systems, such an analysis and simulation 
tool could significantly increase the effectiveness, ef- 
ficiency and reduce the design and development cycle 
associated with entry missions. 
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Fig. 1 Ground tracks of all 6 entry flight trajec- 
tories 



Fig. 2 Comparison of altitude and velocity for 
Mission A1 
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Fig. 3 Comparison of bank angle and a for Mis- 
sion A1 Fig. 6 Comparison of bank angles with different 

nominal a-profiles for Mission A1 



Fig. 4 Comparison of altitude and velocity for 
Mission A2 and A3 
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Fig. 7 Comparison of trajectories with different 
nominal a-profiles for Mission A1 




Fig. 5 Comparison of bank angle and a for Mis- Fig. 8 Heat rates with different nominal a- 
sion A2 and A3 profiles for Mission A1 












